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Abstract 

The distribution is calculated of the electron spin polarization under 
current-driven spin injection from a probe to a ferromagnetic film. It 
is shown that the main parameters determining difference of the spin 
polarization from the equilibrium value are the current density and the 
spin polarization of the probe material, while the relation between the 
probe diameter and the spin diffusion length influences the result very 
weakly, to a certain extent. A possibility is shown of reaching inverse 
population of the spin subbands at distances from the probe boundary 
comparable with the spin diffusion length. 



1 Introduction 

The current-induced spin injection, i.e., appearing nonequilibrium spin 
polarization near boundary between two conductors under current flowing 
through PQ, is one of the main spintronic effects together with the effect of 
spin torque transfer from conduction electrons to lattice [21 [5] . The spin 
injection manifests itself as breaking the thermal equilibrium between spin 
subbands. 

Nowadays, the spin torque transfer effect has been studied with much 
more details than the spin injection. However, an attractive problem oc- 
curs that relates with reaching high injection level. We mean a possibility 
of creating inverse population of the spin subbands in a ferromagnet with 
laser effect in THz and IR ranges (10 12 -10 14 Hz) g]-[9]. The main obsta- 
cle for this idea is necessity of high current density j > 10 9 A/cm 2 . 

A scheme was proposed for reaching high current density in a microprobe- 
thin film system 10 . If the film thickness h is small compared to the probe 
radius R, then the current density in the film near the probe is R/2h times 
the current density in the film. So a problem appears of calculating spin 
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injection in such a system as a function of the system parameters and the 
current value. Such a problem is considered in the present work. 



2 The model and main equations 

Let us consider a ferromagnetic film of h thickness which current / is 
carried to by means of a cylindrical probe of R radius. The current in the 
film is spin-polarized, so that its flowing is accompanied with appearing a 
nonequilibrium spin polarization near the probe at distances comparable 
with the spin diffusion length I. 

The spin polarization is defined as 

P = !H^, (1) 
n 

where n + and n_ are the partial densities of the conduction electrons 
with spin magnetic moments aligned parallel and antiparallel to the crystal 
lattice magnetization vector, respectively, n = n++n- is the total density 
of the conduction density, which is assumed to be constant. It follows from 
Eq. © 

77 

n± = ~(l±P). (2) 

The spin polarization distribution under steady conditions is deter- 
mined by the spin current continuity equation 

VJ = -^^, (3) 

where P is the equilibrium spin polarization, r is the relaxation time of 
the longitudinal (collinear with the lattice magnetization) component of 
the electron spin polarization, 

J = |(j+-j-) W 

is the spin current density 

j± = e/i±n±E - eD±Vn±, (5) 

are the partial current densities created by the electrons with two opposite 
spin directions, E is the electric field, /j,± and D± are the partial mobilities 
and the partial diffusion constants, respectively. 

The electric field E can be expressed via the total current density 
j = j+ + j_. The spin current density Q takes the form 



2e 

where 



H {Q(P)j-enD(P)VP}, (6) 



r>(P\ _ ^+ ~ M- + (/*+ + fi-)P , . 

D(p) = V+D- + + (M+£>- ~ V-D+)P m 
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Substitution of Eq. (|6]) to Eq. (|3} gives an equation for the spin polar- 
ization P which is substantially nonlinear as follows from Eqs. (J2J, (Q 
and Besides the explicit linear- fractional dependence of Q{P) and 
D(P) coefficients on P, the partial mobilities and diffusion constants of a 
degenerate electron gas (metal) depend in general on the Fermi quasilevels 
of the spin subbands and, subsequently, on the partial densities n±, which, 
in their turn, are expressed via spin polarization P (see Eq. l[2)|). The form 
of that dependence is determined by many factors, such as the form of the 
spin subbands (dispersion law), the carrier scattering mechanism, etc. 

In many cases, such as the problem of the switching magnetic con- 
figuration by spin-polarized current, a linear approximation in the cur- 
rent density j and proportional to it nonequilibrium spin polarization 
AP = P — P appears to be sufficient. In this approximation, the spin 
polarization P in Q(P) and D(P) coefficients is replaced with its equi- 
librium value P, so that the coefficients mentioned take constant values 
Q = Q(P) and D = D(P), and Eq. <(3j becomes linear equation, namely, 

V 2 P-^^ = 0, (9) 

where I = V Dr. 

The linear approximation becomes invalid under high spin injection 
corresponding to the spin subband inverse population, when P < 0, so 
that AP < 0, |AP| > P 8 . However, the situation is simplified notice- 
ably if it is supposed that the carriers in both spin subbands have the 
same mobilities and diffusion constants, /i_ = fi + = fi, P_ = D + = D. 
In such a case, we have Q(P) = P, D(P) = D. The substitution of Eq. ([6]) 
into Eq. © with the electric charge conservation condition Vj = taking 
into account gives the equation 

v 2 p-j^_£_£ = , (io) 

where 

I = \/~Dt % j D = enD/l = enl/r. (11) 

The form of this equation depends on neither the carrier degeneration, nor 
the dispersion law, nor the scattering mechanism. The solution of such 
simplified problem, without having any pretension to obtaining quanti- 
tative results for particular materials, allows to find a qualitative picture 
and estimate the orders of magnitude. 



3 Spin polarization distribution 

Under spin injection conditions, the spin polarization differs from its equi- 
librium value at the distances from the injector comparable with the spin 
diffusion length I. If the film lateral size is large in comparison with that 
length, the current density distribution may be considered as axially sym- 
metrical one irrespective of the geometry of the other electrode closing the 
electric circuit. It follows from the electric charge conservation condition 

Vj = ™(rj)=0 (12) 
r ar 
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that the current density distribution in the film near the probe takes the 
form 

^=^7 = 2^' (13) 

where r is the distance from the probe axis, the other notations being 
indicated above. 

The substitution of Eq. (|13|l into Eq. (|10|l gives the following equation 
in polar coordinates: 

d 2 P ,„ „ ,ldP P~P „ 

where 

■Hit*. < 15 > 

The complete solution of Eq. ((14} has the form [11] 

^) = My)>My)+^My)]- ( 16 ) 

where /„, if„ are the modified Bessel functions of the first and second 
kind, respectively, Ci, C2 are the integration constants. 

It follows from P(oo) = P boundary condition, that C\ = 0. The C2 
constant can be found from the spin current continuity condition at the 
boundary between the probe and the film 

J(R) = l^j(R) P (R)- jD l^ 

= ^Qij(R) (Mi • M(i?)) , (17) 

where Q = — ^— is the the probe spin polarization (a± being the 

u+ + a- 

partial conductivities in the probe), M(7?) is the unit vector along the 
film magnetization near the probe, Mi is the same quantity for the 
probe. To obtain inverse population of the spin subbands in the film, 
Mi ■ M(i?)^ < condition is necessary. 
The substitution of Eq. (JTBJ into Eq. gives 

a (ft,. ft^-PjM (18) 

so that the electron spin polarization in the film near the probe takes the 
form 



P(r) = P + {Qi (Mi ■ M(R)) -P} J ^ (^Y KV ^}\ - ( 19 ) 



It follows from Eq. (|19[1 that the spin polarization approaches monotonously 
to the equilibrium value with increasing the distance from the probe 
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Figure 1: Spatial spin polarization distribution near the probe at R/l = 
20, j(R)/jD = 1 and various values of Qi/P ratio. 



(Fig. [TJ. The maximal negative value of the nonequilibrium spin po- 
larization AP is reached at the probe boundary, 

K "(f) 

AP(R) = {q 1 (m 1 -M(R))-P\ J -^ (20) 

As to the dependence of the nonequilibrium spin polarization on the 
current density, it is necessary to have in mind that the current den- 
sity j appears in Eqs. (|19[) and (|20[) not only as an explicit factor, but 
also in v parameter (see definition JTSJl). Because of such a reason, 
the spin polarization near the probe tends to a limiting value, P(R) — * 

Qi -M(P)^. Note, that the nonequilibrium spin polarization ap- 
pears in the case of a nonmagnetic probe also (Qi = 0), that means re- 
placing the spin-polarized electrons with non-polarized ones in the vicinity 
of the probe. In that case, the spin polarization remains positive and tends 
to zero at high current density (AP — > — P, P — > 0). 

In Fig. [T] a coordinate dependence is shown of the spin polarization 
(referred to its equilibrium value) at given values of j(R)/jD and R/l 
ratios and various values of Qi/P ratio. The spin polarization near the 
probe becomes negative at high enough values of the latter ratio. 

In Fig. [2] the spin polarization on the probe boundary is shown as 
a function of the current density at given value of R/l ratio and various 
values of Qi/P ratio. 
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Figure 2: Spin polarization at the boundary between the probe and the film 
as a function of the (dimensionless) current density j(R)/jD at R/l = 20 and 
various values of Qi/P ratio. 



As numerical analysis shows, the results weakly depend on R/l ra- 
tio. The parameters which influence the spin injection substantially are 
j{R)/jD and Qi/P ratios. 

At n ~ 10 22 cm" 3 , I ~ 1CT 6 cm, r ~ 1CT 12 s, we have j D ~ 10 9 
A/cm 2 . With the current density in the probe jo ~ 10 s A/cm 2 (such a 
value was reached, e.g., in [13]) and R/2h > 10 (this ratio was equal to 5 
in [12]), we find j(i?) ~ 10 9 A/cm 2 ~ ju, that is sufficient, in accordance 
with Fig. [2] for negative spin polarization near the probe at Qi/P > 1. 

4 Current-induced magnetic field and sd 
exchange effective field effects 

At high current densities, influence becomes substantial of the current- 
induced magnetic field (the Ampere field). For rough estimates, we may 
use a simple formula for the magnetic field of a straight current 

H=L (21) 

without factor 2 of current ("the half of a long wire"). The substitution 
of I = 2n Rhj(R) into Eq. (|21jl gives the magnetic field at the probe 
circumference with given current density: 

H = 2irhj{R). (22) 
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At h = 10 nm and j(R) = 10 9 A/cm 2 we have H « 600 Oe. 

If the film anisotropy field H a is lower than the latter value, then the 
lattice magnetization near the probe is directed along the probe circum- 
ference. To align the probe magnetization opposite to the film magneti- 
zation, the probe anisotropy field H a i is to be larger than the indicated 
H value, and the probe is to be magnetized by such a current of opposite 
direction that induces H > H a i field at the probe circumference. 

The inverse population is prevented also with the sd exchange interac- 
tion of the injected electrons with the magnetic lattice of the film, which 
tends to align the electron spins in the film parallel to ones in the injector. 

The sd exchange energy is 



U sd = 



t J (M(r) ■ m(r)) d 3 r, (23) 



where a is the dimensionless constant of the sd exchange interaction, 
M = MM is the lattice magnetization, 

m = fi B nP(v)iA (24) 

is the electron magnetization in the film, [ib is the Bohr magneton. 

Substitution of Eqs. ((24} and (fl9|) into Eq. (J23J gives the following- 
expression for the nonequilibrium part of the sd exchange interaction en- 
ergy: 

OO 

AU sd = -afi B nM ■ 2nh J {P(r) - P} r dr = 

R 

= -a^griM {Qi (Mi • M) — P\ ■ 2-kRM. (25) 

3d 

The corresponding sd exchange effective field 

has an order of magnitude 

H sd ~/x B anQi^. (27) 

JD 

This field becomes comparable with the molecular field ~ 10 6 Oe at 
j(R) ~ Jd- To prevent switching antiparallel initial configuration ( (Mi ■ M(R) 
0) to parallel one, it is necessary to pin the magnetization of the film by 
means of induced anisotropy with the aid of an antiferromagnetic sub- 
strate. 



5 Conclusion 

The analysis within a simplified model shows general possibility of reach- 
ing nonequilibrium negative spin polarization (the spin subband inversion) 
with using probe/film structures. However, realization of this possibility 
needs some action to prevent switching the antiparallel configuration to 
parallel one. 
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